ANNALES 
UNIVERSITATIS MARIAE CURIE-SKLODOWSKA 
LUBLIN—POLONIA 
VOL. XXXIX/XL, 10 .' SECTIO AA n 1984/1985 


Zakład Chemii Teoretycznej 
Instytut Chemii UMCS 


Andrzej DĄBROWSKI 


Równanie izotermy adsorpcji z dwuskładnikowych ciekłych mieszanin na ciałach 
stałych uwzględniające energetyczną heterogeniczność adsorbentu oraz rozmiary czą- 
stek składników roztworu i jego numeryczna weryfikacja 


An Isotherm Equation for Adsorption from Binary Liquid Mixtures on Solids Involving Sur- 
face Heterogeneity and Differences in Molecular Sizes of Components and its Numerical 
Verification 


- VpasHenue n30TepMbr a4COpÓLMU M3 64HapHbIX. 2KM7IKIEX. CMeceit Ha IOBEPXHOĆTHM TBep/IbIX TEN 
yuwrWBalou(Me 3HepreTH"ecKyio reTeporeHHOCTb HOBepXHOCTM.  pa3Mepbl MO.IEKYJI M ero Hy- 
MepuieckaA mnpoBepka 


RZ i i WPROWADZE NIE 


w ostatnich latach opublikowano szereg prac dotyczących wpływu 
powierzchniowej heterogeniczności na proces adsorpcji z roztworów 
na ciałach stałych [1-12], W monografii [13Jzamieszczono najważ- 
niejsze osiągnięcia z zakresu tej tematyki, z uwzględnieniem zarówno 
mono» jak i wielowarstwowych modeli fazy powierzchniowej, Równo- 
cześnie praca ta stanowi zwięzły przegląd najważniejszych osiągnięć 
' dotyczących teoretycznych i praktycznych aspektów omawianego za” 
gadnienia, Z pracy [13] wynika, że stan teoretycznych badań dotycza- 
cych procesu adsorpcji na powierzchniach ciał stałych z roztworów 
składających się z cząsteczek o różnych rozmiarach należy uznać 


za wysoce niezadowalający, Z drugiej strony, większość realnych  . 
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ukiadów adsorpcyjnych wiacza zarówno adsorbenty o. powierzchniach 
heterogenicznych jak i wspomniane wyżej rodzaje roztworów, 

w pracy E141, w ramach formalizmu kinetycznego, wyprowadzo- 
no równanie izotermy adsorpcji monowarstwowej, w którym uw zględnio= 
no zarówno różnicę w rozmiarach adsorbujących się cząsteczek jak 
i enercetyczną heterogeniczność adsorbentu, Tamże wykonano bada= 
A bdidejowe, pokazujące wpływ wow czynników oraz oddziaływań Bo 
między cząsteczkami roztworów na kształt izoterm nadmiarowej adsorpcji 

' j odchylenia układu modelowego od zachowania idealnego, 

Ww ostatniej pracy Rudzińskiego i współpracowników [15], 
w której wykorzystano metodę przybliżenia kondensacyjnego zastosowaną 
po raz pierwszy w teorii adsorpcji z roztworów na heterogenicznych 
powierzchniach ciał statych przez Ościka i współpra= 
cowników C16, 173 oraz Dąbrowskiego i wspótpra- 
cowników [2], wyprowadzono analityczne równanie, którego postać 
jest zbliżona do wyrażenia zaproponowanego w pracy C141. 

W przeciwieństwie do równania wyprowadzonego w pracy L14J 
równanie Rudzińskiego nie redukuje sie do ogólnie znanego równania 
Everetta, opisującego proces adsorpcji z roztworów na homogenicznych 
powierzchniach ciał stałych L181. Wynika to niewątpliwie z faktu, że 
metoda przybliżenia kondensacyjnego daje zadowalające rezultaty w. 
odniesieniu do procesu. adsorpcji na bardzo silnie heterogenicznych 
adsorbentach [ 2-121, Należy przypuszczać, że w tym ostatnim przy- 
padku wyrażenie zaproponowane przez Rudzińskiego i współpracowni= 
ków może być bardziej użyteczne. Jednakże celem prezentowanej pra- 
cy nie jest porównanie obu wspomnianych wyrażeń, Porównanie takie 
będzie zaprezentowane w następnej publikacji z tego cyklu C191. 

Podstawowym celem niniejszej pracy jest bowiem zbadanie 
właściwości równania wyprowadzonego w pracy £141i potwierdzenie 
jego użyteczności dla opisu ciekłej adsorpcji na słabo i silnie hete- 
rogenicznych powierzchniach ciał stałych. Realizacja tego celu będzie 
osiągnięta poprzez przeprowadzenie opisanych nizej symulacji nume- 
rycznych. Dła prostoty rozważań wprowadzimy model IBP układu 

| adsorpcyjnego (ideal behaviour in both adsorbent- and bulk phases) 
[14], Zgodnie z wcześniejszymi rozważaniami [. 141, w takim przypad= 
"ku rodzaj topografii miejsc adsorpcyjnych nie ma znaczenia, Wobec 

: faktu rozważania adsorpcji z roztworów złóżonych z cząsteczek o 
różnych rozmiarach $tąsowane będą. ułamki objętościowe w miejsce. 
zwykle wprowadzonych ułamków molowych, 
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ROZWAŻANIA OGÓLNE 


Jedno z podstawowych równań opisujących proces adsorpcji z 
roztworów na powierzchniach heterogenicznych, charakteryzowanych 
przez ciągły rozkład energii adsorpcji, ma następującą postać L 131: 


4 - | 93 (2,01) P eeu a) 
R 


gdzie E jest różnicą potencjałów adsorpcyjnych obu składników 


12 
roztworu , $i oznacza średni ułamek objętościowy składnika "1", 
Qi dotyczy adsorpcji.na homogenicznym płacie powierzchni charak= 

teryzowanym przez różnicę energii adsorpcji E42 natomiast S2 jest 


przedziałem możliwych zmian wiełkości E4 $i jest ułamkiem obję= 


, 
tościowym składnika "1" w fazie Fco ti t a Równanie (1) rozwia- 
zuje się zwykle metodą transformat Stieltjesa C 131. Metoda ta nie 
może być jednak stosowana w przypadku, gdy roztwór Składa się z 
cząsteczek o różnych powierzchniach siadania, Wtedy bowiem izoterma 
lokalna 9 > ma postać funkcji uwikłanej L 1831. 

Załóżmy teraz, zgodnie z rozważaniami pracy L141, że hetero- 
geniczna powierzchnia ciała stałego składa się z M miejsc adsorpcyj- 
nych rozmieszczonych na N homogenicznych płatach, Spetniony jest 
przy tym następujący warunek: | 


N 
M a i2. My « 
dem rs 


aż! 


gdzie M. jest liczba miejsc adsorpcyjnych k-tego typu, Wówczas 


zamiast równania (1) możemy napisać następujące wyrażenie L 141: 


s A. s. 

$4 = 26 fc 9 1k vg l (2) 
/ ks 

' gdzie fi =.M,/M  oznacza ulamek miejsc adsorpcyjnych umieszczo= 

nych na k=tym płacie, związanych z energią E12,k* Odpowiednia stała 

równowagi adsorpcyjnej ma. następującą postać C14, 187: 
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v (4t) SK jk | (3) 


gdzie ©. = 1 -01. 4 2T 1% d * Ponadto 


Ke = exp (E45 I RT) ` (4) 
| oraz ; 
E "«rE,,=E ; i (5) 


12,k 1,k 2,k 


Parametr Ek (i=1,2) oznacza energię adsorpcji składnika "i" na 
k-tym płacie. Symbol r występujący w równaniu (3) jest stosunkiem 
pól powierzchni siadania składników 2 i 1, t r = wawi . Na pod- 
stawie równania (2), przy użyciu metody kinetycznej opisanej w 
pracy [14], wyprowadzono następujące równanie (6): 


(D * 


EK = exp (E, IRT) - "E "(60)" (6) 
Ą 2 : 1 


gdzie Ki Ei. oznaczają odpowiednio średnie stałych K, oraz energii 


E12 k i dotyczą całej heterogenicznej powierzchni ciała stałego [14]. 
* 


Parametr c będący parametrem heterogeniczności można określić w 


sposób następujący [5]: 


c= 1/(e, + 1) . c (7) 


gdzie €4»9 jest parametrem wolnym, W konsekwenc ji c«( 0,1 > 
i charakteryzuje kształt quasi-gaussowskiej funkcji rozkładu energii 


- adsorpcji, Gdy c=1, równanie (6) redukuje się do postaci równania 


Everetta opisującego proces adsorpcji na homogenicznej powierzchni 
ciała stałego z roztworów złożonych z cząsteczek o różnych roz= 
miarach C 181, Gdy r=1, to równanie (6) można przekształcić, otrzy” 
mując dokładne analityczne rozwiązanie wyrażenie (1)C[ 21. 
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ROZWAZANIA TEORETYCZNE 


Przyjmijmy następujący, typowy rozkład Gaussa zmiennej 
E LE rE „-E s 


12 2 
3 (E P E? )? 


gdzie A jest parametrem charakteryzującym szerokość wprowadzone= 
go rozkładu, natomiast EJ „odpowiada wielkości E„„, dla której funkcja 


F(E,2) osiąga wartość: maksymalną, Cdy parametr A dąży do 0, to 
wtedy funkcja F(E.2) przyjmuje postać funkcii $ -Diraca, Funkcja 
ta charakteryzuje omawiany proces na homocenicznym płacie powierz- 


chni związanym z energią adsorpcji E = EQ . Funkcja F(E,,) 


12 "12 
określona równaniem (8) spełnia następujący warunek normalizacji: 


4 o 


ne da) A Brae Lot (9) 


— 


o 
12 
znormalizowane funkcje F(E,,). Dla celów porównawczych wszystkie 


Ryc, 1, przedstawia obliczone dla E,, = 0 i różnych wielkości A 
pokazane tu funkcje zostały narysowane w jednakowej skali, Z oma- 
wianego rysunku wynika, że dystrybucja | F(E,,) związana z parame- 
trem A =a 0,25 (krzywa a) przybliża funkcję © -Diraca i odpowiada 
12 < 70:65 kcal, + 0,65kcal^. 


Innymi słowy, omawiana funkcja charakteryzuje powierzchnię adsorben- 


jej nieznaczny przedział zmian energii: E 


tu bliską homogenicznej, Funkcja F(&4,) związana z parametrem 
Aa (krzywa d) koresponduje z szerokim zakresem zmian energii: 
E12 < =5.5kcal,+5,5kcał > i charakteryzuje wysoce heterogeniczna 
powierzchnię ciała stałego, 

Uwzględniając wyrażenie (2) należy zastąpić ciągły rozkład (8) 
rozkładem dyskretnym, W jednej z ostatnich prac Czarniecki i Jaro= 
niec £ 20 J wykazali, że taka procedura jest w pełni uzasadniona przy 


uwzględnieniu dostatecznie szerokiego przedziału zmian energii E12* 
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Ryc, 1, Znormalizowane funkcje rozkładu określone równaniem (8) 
obliczone dla różnych wielkości ^ oraz E: = 0: A = 0.25 (a); 
Am 0.5 (b); A = 1 (c); A -2 (qd) 
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W przypadku dyskretnego rozkładu energii adsorpcji warunek normali- 


zacyjny (9) może być przepisany w następujący sposób: 


N 
; CY | (10) 
| k»1 
gdzie 
Wei 
& = P (Eiai) / 0. F (E12) (11) 
i EAT 


D 


Należy podkreślić, ze znaczenie parametru í , wynikające z równania 


(11) jest równoważne definicji (2). CZE cel naszych dal- 
szych rozważań, na ryc, 2 przedstawiono funkcję P(E 12! obliczona 
przy założeniu, ze Â « 1 oraz E15" 0. W tym przypadku przedział 
znaczących zmian wielkości E12 jest następujący: E 426 € -3.2kcal,+ 
+3,2kcal > . Przedział ten: został podzielony na 16 równych co do 
wielkości podprzedziałów, co oznacza, że wartość górnego indeksu 
w sumowaniu . (10) jest równa: N=17, W celu wykazania użyteczności 
równania (6) zostały wykonane opisane niżej operacje numeryczne: 


i 1, Przy użyciu wzoru (12) obliczano ułamki objętościowe et 
i $ neus 
1 
1 
+ r(1-xl) 


4 - 
4 a 
s 
dla danych wielkości r oraz xy» 
2, Stosując równanie (4) obliczano stałe równowasi K, Opisują- 
ce adsorpcję na homogenicznych płatach związanych z-energiami 
E 42,k oszacowanymi dla danej funkcji. F(E42) (patrz takze "ye 2). 
3, Stosując równanie (2) obliczano ułamki objętościowe oi Š 
-W tym celu na podstawie równania (11) znajdowano ułamek miejsc 
adsorpcyjnych umieszczonych na k=tym IAW OWIÓSY płacie, charak= 


ZE przez energie adsorpcji „ Ułamek objętościowy 


E12, k 
Qi 1,k odpowiadający adsorpcji na tym samym płacie obliczano roz= 
wiązując numerycznie równanie (3), w którym występujące wielkości 


stałych Ki znajiowano na podstawie równania (4). 
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Ryc. 2, Schemat ilustrujący numeryczną procedurę stosowaną w celu 
weryfikacji równania (6) 


4. Przeliczano ułamki objętościowe $1 na ułamki molowe xi. 
stosując następującą formułę: : 


"AGS 


The r$ + (1-95) (13) 


a następnie znajdowano izotermę nadmiarowej adsorpcji Nan (x1), 
opisującą omawiany proces na całej heterogenicznej powierzchni ciała 


stałego, Stosowano przy tym następujący wzór: 


e 1 
X SĄ. - 3E 

Nf. z z 1 (14) 
r+X; (1-r) 


Symulowane numerycznie izotermy $4 = CZE oraz Ne NĄ (x) 


traktowano jako izotermy doświadczalne, 


5. Stosując liniową postać wyrażenia (6) BO vx eks- 
perymentainą funkcję pi = CACORE 


` 
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s Ay 1.9! 
in ail - € in ( £)71 + c In aż 


ECUN qu 


A (15) 


6. Następnie, według równań (6),(12) - (15) obliczano teoretycz- 
ne izotermy nadmiarowej adsorpcji. Aby rozwiązać równanie (6) wpro” .: 


wadzano parametry c oraz K znalezione na podstawie równania (15), 


7, Obliczano sumę kwadratów odchyleń (5KO) zdefiniowaną w 
następujący sposób: i 


1 7 
2 
SKO = X NyeKSP_ u weteon | (16) 
i=1 | 


gdzie I oznacza całkowitą liczbę punktów eksperymentalnych, Symbole 
"SE Pia oraz BZ: oznaczają doświadczalne i teoretyczne 
wielkości izoterm nadmiarowej adsorpcji, Wykonując obliczenia nurae- 
ryczne według punktów 1-7 parametr r w odpowiednich wzorach zapi- 
sywano w postaci: p=ljr, » gdzie r.= wą lwa . W trakcie obliczeń 
przyjęto, że F=0,6, 0,7, m, 1.4, xh=0.05, 0,15, „w 0,95, Is10, N=17, 
RT=0.6kcal/mol. 


WYNIKI OBLICZEŃ NUMERYCZNYCH I ICH DYSKUSJA 


, Na ryc. 3 przedstawiono liniowe zależności ; 

In E: -AS d vs. in E (1-$))/($4) 1 à obliczone zgodnie z 
równaniem (15) dla E170 oraz Â 20,25 (prosta c) i A =1 
(prosta a oraz b), Dła obu wielkości A punkty eksperymentalne 
znajdowano stosując równanie (2) dla r40.6 (białe kółka) oraz 
Vj=1.4 (czarne kółka), Z ryc. 3 wynika, że równanie (15) będące 
liniową postacią równania (6) dobrze reprezentuje symulowane nume= 
rycznie dane adsorpcyjne zarówno dla A = 0.25 (mała powierzchnio- 
wa heterogeniczność ) jak i dia A -1 (znaczna powierzchniowa 
heterogeniczność ) (patrz także ryc, 1). W pierwszym przypadku pa= 
rametr c określany bezpośrednio przez nachylenie linii prostych jesł . 


praktycznie niezależny od różnicy molekularnych rozmiarów obu skład 
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Ryc. 3. Liniowa zależność obliczona według równania (15) dla symu- 


"Jowanych danych doświadczałnych, E1570. Pozostate parametry ad- 


sorpcyjne: A 21; F17*1.4 (ciągła prosta a i czarne kółka), A =1; riti 


(prosta c i czarne kółka), Â =0.25; r,-0.6 (prosta c i białe kółka) 


1 


ników roztworu (c=0,997 dla r,*9.6 oraz c=0,998 dla r4=1.4), lecz w 
ostatnim przypadku omawiany parametr c jest nieznacznie zależny od 
r4 {c=0,643 dla ra 


padku adsorpcji na powierzchniach heterogenicznych z roztworów zło= 


"0.6 oraz ce0,?75 dla r,-1.1), Innymi słowy, w przy= 


żonych z cząsteczek o różnych rozmiarach parametr ra zaburza praw= 
dziwą wartość paremetru heterogeruczności c. Na pierwszy rzut oka 
wniosek taki wydaje się zaskakujący, ponieważ bez względu na zmiany 
parametru ri parametr c charakteryzujący za pośrednictwem parametru 
ŻA «const kształt funkcji rozkładu (8) powinien mieć wielkość stałą, 
Niemniej powyższy wynik pozostaje w pełnej zgodności z badaniami 


modelowymi wykonanymi na podstawie równania (6) i prezentowanymi 
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w pracy [ 147] (patrz ryc, 1 i praca [ 1427). W cytowanym artykule 
stwierdzono co następuje: "s, stosunek pól powierzchni siadania zaje 
mowanych przez oba składniki może wzmagać się dla r < 1 lub ob- 
niżać. dla r4 > 1 wpływ powierzchniowej heterogeniczności”, W obli- 
czu powyższego spostrzeżenia należy, stosując równanie (6) do opi= 
su konkretnych danych doświadczalnych, wprowadzać stosowane po- 
prawki ze względu na parametr c (patrz dalsze rozważania). Z dru= 
giej strony, gdy adsorpcja zachodzi na słabo heterogenicznych po- 
' wierzchniach ciat statych, wtedy wplyw parametru r, na c moze byc 
zaniedbany, Jest rzeczą interesującą, że stała równowagi adsorpcyj- 
nej.K, podobnie jak parametr c, jest praktycznie niezależna od EU 
gdy c dąży do jedności, Na przykład ^ 40,25, ce1, K=1,008 dla 
250,6, i K=0,996 gdy r121.4. Jednakże różnice w wielkościach K 
wzrastają równocześnie z EU gdy parametr A rośnie, 

Powyższe wnioski są w pełnej zgodności z wynikami przedsta 
- wionymi na ryc, 4. W części B tego rysunku przedstawiono liniową 
zależność (15), którą zasiosowano do opisu symulowanych numerycznie, 
na podstawie. równania (2) dla ŻA =2, punktów eksperymentalnych, 
W części A ryc. 4 przedstawiono odpowiednie eksperymentalne i teo- 
retyczna izotermy nadmiarowej adsorpcji, Oznaczenia na ryc, 4 są 
takie, jak na ryc, 3, z tym że linia przerywana z przekreślonymi 
kółkami (ryc. 4) odpowiada sytuacji, gdy r,"1. Parametry c, R oraz 
SKO odpowiadające wynikom przedstawionym na ryc, 3, 4 są zamiesz- 
czone w tab, 1. Z tabeli tej wynika, że bez względu na zmiany para- 
metru Â wielkości K dla r,"1 odpowiadają energii E1570» przy czym 
w tym przypadku funkcja (8) osiąga wartość maksymalną, tj. Kal. 
Przedstawiona wyżej uwagę można wytłumaczyć jak następuje: gdy 
rtl to równanie (6) może być przekształcone do następującej pos- 
taci: 


In K -LE - Ej) Rri « ET IRT (17) 


Ponieważ funkcja (8) jest symetryczna ze wzgiędu na swój argument, 
wielkość Eiz jest wielkością średnią różnicy potencjałów adsorpcyj- 
nych obu składników, dla których funkcja Gaussa osiąga wartość 
maksymalną, tj, ET, "E1570 ikea, Ponadto wynik ten dowodzi do- 
kładności procedury polegającej na zastąpieniu całki (1) sumą (2). 

` Zauważmy, że równanie (6) dla c=l oraz r,=1 przyjmuje posłać 


najprostszego równania Everetta [221], które dla ŘKR=K=1 przewiduje 
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wt) 


in (1-9 


0 E 
in [0 - 652/091] 


Ryc. Symulowane izotermy nadmiarowej adsorpcji (A) oraz odpo- 


wiadajace im funkcje (15) (B) obliczone dla ES 


 Kólkami oznaczono punkty eksperymentalne, natomiast linie oznacza- 


279 oraz =2, 


ją funkcje teoretyczne. Pozostałe parametry adsorpcyjne mają wartoé- 
ci; r4*1,4; (linia ciągła a i czarne kółka): r;71 (linia przerywana b 
i przekreślone kółka); F170.6 (linia przerywana c i białe kółka) 


zerową wielkość izotermy nadmiarowej adsorpcji w całym przedziale 
stężeń, Jak wynika z rye. 4 oraz z pracL 2-111 w trakcie adsorp- 
cji z roztworów na heterogenicznych powierzchniach ciał stałych izo- 
' ferma nadmiarowej adsorpcji jest rózna od zera, nawet wtedy, gdy 
Rel oraz ryth Analogiczny wniosek jest prawdziwy w odniesieniu 
do adsorpcji na powierzchniach homogenicznych, gdy c=i oraz r,t 1 
£191. : 

Sumujac uzyskane rezultaty możemy stwierdzić, że testowane 
równanie (6) zapewnia dobrą reprezentację danych doświadczalnych 
w szerokim zakresie zmian powierzchniowej heterogeniczności ciała 
stałego, przy czym parametr c uzyskany na podstawie tej formuły za- 
leży od Fa gdy parametr A rośnie, Prawdziwą wielkość parametru c 
można znaleźć za pomocą diagramu ZA = A (c) podobnego do tego, 
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Tab, 1, Parametry adsorpcyjne równania (6) obliczone dla E) 20 
i 12 
oraz A^ a 0,25, 1, 2 
DĄ: si > NE t x T 
SA r4” (r) I K SKO Nr ryc. 
0.6 0.978 ^^ 1,009 0,00009 
0.25 1,0 0,983 1.001. ^ ^ 0,000006 3 
14 :' . 0,988 0,998 0,00004 
0.6 0,643 1.141 | 0,0002 
1,0 | 0,679 | 1,001 0,0000 1, 
1.5.5 2:4 (| 9??75.. 110.935 0,0004 > 
0.6 0.286 1.821 0,0004 
2 1.0 0,345 1.000 0,00003 4 


14 0.422 0,715 0,0005 


" 
i 


w 


— 


0 E : 
N 0 0.2 0.4 0.6 0.8 E 10 


A 
Ryc 5, Zależność ^ PAS (c) obliczona przy użyciu równań (2) oraz 
(15) dia $0. Pozostałe parametry adsorpcyjne mają następujące 
wartości: rQ71.4 (linia ciągła i czarne kółka); r.=1 (linia przerywana 
"i przekreślone kółka); P1=0,6 (linia przerywana i białe kółka) 
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jaki zaprezentowano na ryc, 5. Rysunek ten potwierdza nasze uprzed- 
nie wnioski dotyczące wzoru (6), a mianowicie 147: 

1) równanie (6) koresponduje z gaussowską funkcją zmiennej 
E 12* A 

2) równanie (6) przewiduje, że parametr heteroseniczności 
cec (Ô ) należy do przedziału: c € (0.1), 

3) edy eml, to parametr Â =O i funkcja (8) przyjmuje postać 
funkcji $ -Diraca, 

Jednocześnie, jak wynikaz z ryc. 5, równanie (6) daje nieza- 
dowalające wyniki, gdy Zw >25 i c 0.3, ponieważ wówczas wykres 
A =»A(c) jest mało zależny od zmian A . Jednakże parametry c, 
cliarakteryzujące większość układów eksperymentalnych są większe 
niż c=0,5 £3, 231. l 

Aby obliczyć funkcje (8) charakteryzujące dowolný układ 
eksperymentalny wymagana jest znajomość parametru ES e Procedura 
służąca temu celowi oraz zastosowanie równania (6) do analizy 
rzeczywistych układów adsorpcyjnych, a także inne zagadnienia do- 
tyczące ciekłej adsorpcji z roztworów złożonych z cząsteczek o 
różnych powierzchniach siadania będą tematem następnej publikacji 
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SUMMARY. 


This paper presents the numerical verification of equation des 
cribing an adsorpłion process from solutions composed of molecules 


of different sizes on. heterogeneous surfaces of solids, 


B sro cTaTbu npeAcTaBJeHa HyuepurecXaq nponepka ypaBHe- 
HMA OIMCHBARDUero.tnüpouécc 8JCcOpÓ:4M MƏ pacTBOpoB COCTOALIUX M3 
MOJeKyX paaHoro paamepa Ha reTeporeHHHX NoBepXHOCTAX T3EpNAX 
Tea. 
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